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Abstract 

In this paper we consider two statistical hypotheses for the families of Wishart type 
distributions. These distributions are analogs of the Wishart distributions defined 
and parametrized over a Lorentz cone. We test these hypotheses by means of maximal 
invariant statistics which are explicitly derived in the paper. The testing problems, 
respectively, concern the hypothesis that parameters are in a sub-Lorentz-cone, and 
the the hypothesis that two observations have the same parameter. 
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1. Introduction 

The primary aim of this paper is a detailed study of maximal invariant statis- 
tics for a family of Lorentz type Wishart distributions, or simply a Lorentz Wishart 
model. A Lorentz type Wishart distribution, similar to the classical Wishart distri- 
bution [11], naturally arises as the distribution of the empirical sample covariance 
of a multivariate normal distribution [8], or directly as the joint spectral density 
of observable variables [7]. Generally, by employing the theory of Euclidean simple 
algebras, the Wishart distributions can be defined on each irreducible symmetric 
cone [5], or more generally, on each homogeneous cone, by employing the theory of 
Vinberg algebras [2]. According to the classification of irreducible (symmetric) cones 
[6] there are five types of irreducible cones. The first there types of irreducible cones 
are, respectively, PD„(]R), PD„(C) and PD„(EI), i.e., the cone of positive definite 
matrices over the field of real, complex and quaternion numbers. The forth type of 
irreducible cone is the Lorentz cone (or Mikowski cone), and the fifth type is the 
exceptional cone over the Octonion O. Therefore, correspondingly, the three first 
types of Wishart distributions are real, complex and quaternion Wishart distribu- 
tions which have been well studied in statistics and probability literature. In [8] 
Jensen considered statistical hypotheses for Lorentz Wishart models and obtain a 
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complete solution to the problems of maximum likelihood inference. 



In this paper. ^ we explicitly derive maximal invariant statistics for testing two in- 
variant statistical hypotheses for the family of Lorentz type Wishart distributions. 
The first hypothesis is studied in [8] too, but our approach is more general and it 
highhghts the role of the maximal invariant statistic. The second hypothesis is an 
analog of Bartlett's test. 

The organization of this paper is as follows. In §2 we give a precise definition of a 
Lorentz cone £, the description of the automorphism group G of the Lorentz cone 
£, and the definition of the Lorentz type Wishart distribution W^^^. In §3 we test 
the hypothesis that the scale parameter cr is in a Lorentz subcone Co C C To this 
end, first we identify a subgroup Go d G that acts transitively on Co and derive a 
maximal invariant statistic associated with the hypothesis under the action of Go- In 
§4 we test the hypothesis that two observed Lorentz type Wishart distributions have 
the same scale parameter a. For this testing problem we derive a maximal invariant 
statistic associated with the hypothesis under the action of G. 

2. Preliminciries 

2.1. The Lorentz cone 

Let W 7^ {0} be a Euchdean vector space and let 



denote the inner product on W. Set as usual ||w|p := (w,w), for each w G W. 
Consider the symmetric form ^ on R x given by ^'((A, w), (A, w)) := A^ — 
This is a non-singular symmetric from with signature (1, Dim{W)). The symmetric 
form ^ partitions according to the decomposition M x 1^ as 



where V'o e The set 

C := {(A,w) eRxW\ X> 0,^((A,w)) > 0} 
is a symmetric cone, called the Lorentz cone (generated by W). 

^This paper is a revised version of the second part of the author's doctoral thesis [4]. In the first 
part, using a uniform approach, the maximal invariant statistics over irreducible cones are studied, 
^we use the notations R++ and respectively, for the set of non- negative and positive real 



^ty(w,w) = (w,w'), W,w' e W, 




numbers 
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2.2. The automorphism group 

Define for any symmetric form ^ on a Euclidean space V the orthogonal group 

0(*) := {A e GL{V)\ ^o{AxA)^^}, 

the special orthogonal group 

S0(*) {A e 0(1^)1 det{A) = 1} . 

The connected component of the identity in ©(^f), denoted by SO^(^), is a subgroup 
of SO(^). If ^ is positive definite, and hence an inner product on V, then SO^(^) = 
S0(*). 

Let Gw — G denote the connected component of the automorphism group Kut{C). 
Then for the symmetric form ^ defined in §2.1 we have 

where 

with aoiy : ^ M, Slwo : M ^ and Aw :W^W being linear mappings. Note 
that M++ X SO^ acts transitively on £ and for (a, A) e M++ x SO^(M x VT) we have 

(a, A)(A, w) = ayl 

2.3. The Lorentz type Wishart distributions 

Fix an element e G VT. Let ll^i be the orthogonal complement of Me, i.e., Me^. Thus 
= Me © H^i can be identified with M x Wi. Under this identification, the Lorentz 
cone £ is isomorphic to the homogenous cone 

r2{W) := {S = ( ^A: Ai,A2eM++, wieW^i,det(5):=AiA2-||wi|p>0}, 

\Wi A2 J 

studied in [2]. The isomorphism is given by 

/ Ai wi A .A1 + A2 .A1-A2 , , 

Ui A2 J ^ 
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where the inverse mapping is given by 

0-^ : £ ^ r2{W) 
(A.(A.»)) - (^'^'^ (3) 

By using their general theory of the Wishart distributions for homogeneous cones, 
in [2] Andersson et aL derive the Wishart distribution on V2{V) as 

Tfl 

where > — is shape parameter, m — <\im^{W), and E e V2{W) is the expectation 

vector. Therefore, one natural way to define the Wishart distribution on C is via the 
image of yVv2iv){S\r), E) under the mapping 0. 

Definition 2.1. The Lorentz type Wishart distribution, denoted by W^^j, is the image 
of the Wishart distribution Wv2{v)iS\r],(l)~^ (a)) under the mapping (f) in Eq. (2). One 

in — 1 

can check that the Wishart distribution on L with shape parameter rj > — - — and 
expectation a — {X,w) & C is given by 

where fe.z) e Rx on(i*(m,ij) = 2ir°T'r(>))r(>) - =!fi)Tr^'. 
') = f ■ 

Remark 2.1. First note that the density given in Eq. (4) is the same as formula (29) 
in [8]. Moreover, for every irreducible symmetric cone Q, the Wishart distribution 
on Q, denoted by W^^^,, is well defined and given by the density 

dK^i^) = ,1,.^ exp{-ltr(a-^x)}det(x)V7l^(x), (5) 

rn(^) det((7)2'' z 

where r is the rank of ft, and rQ{-) is the gamma function associated with ft (see 
[10] or [3] for detail). Since C is an irreducible symmetric cone of rank 2, Eq. (5) 
directly defines the Wishart distribution on L. 
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3. Testing for scale parameter in a Lorentz subcone 

Suppose that shape parameter r] is known and consider the Lorentz Wishart model 

M^{Wi,,: aeC}. (6) 

The standard theory of exponential families implies that the ML estimator amie of 
cr e £ e exists for any observation {y, z) e £ and is given by 

^mle{{y,z)) ^ {y,z). (7) 

Now suppose that Wq ^ {0} is a subspace of W. Let Cq denote the Lorentz cone 
generated by Wq. Under the inclusion mapping l : jCq ^ C the Lorentz Wishart 
model 

Mo ^ : <Ji e Co} (8) 

is a submodel of M.. Let us consider the hypothesis 

Ho: a e Co vs. H : a e C. (Tl) 

We will test the hypothesis (Tl) by a maximal invariant statistic wc shall derive 
in §3.2. First we identify a subgroup Go d G such that the hypothesis is invariant 
under it. We proceed as follows. 

3.1. The subgroup of G with invariant action on Cq 

Recall the definition of the symmetric form \& and the inner product '^w in §2.1. Let 
be the orthogonal complement of Wq C VF, with respect to We set 

■^w\Wo X Wo, -^wi- — '^w\W^ x W^, and *o := x Wo) x (R x W^o)- 

By these conventions ^ partitions according to the decomposition M x Wo x Wq as 





( V'o 










' 









^ 








Note that and partition with respect to the decompositions IR x Wo and 
W — Wo^ Wq accordingly as 
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Since (/x,wo) G Cq implies that /i^ > ipwoi"^) = i^wi'^), clearly, Cq C M. x Wq C 
R X 1^ is a subcone of £, i.e., jCq C C 

Next we determine the group Gq :— {g e G\g{Co) C Cq}. Note that 

R++ c R++ X SO^(*) 

is obviously contained in Gq. 

Proposition 3.1. Let A e SO^(^'). Then A{Co) C jCq if and only if A partition 
as 

Ao 



A 



with 



(i) AoeOi^o), i.e., *o o (^o x Aq) = *0; 

(ii) A^^ e 0{^w^), i.e., *;^x o (^^x X A^^) = *i^x, 

('mj flo > 0, and 

(iv) (let{Awa) ■ det{A^i.) > 0. 

Proof. Suppose A e SO^(^') acts invariantly on Cq, i.e., A{Co) C Cq. Let 



-1 


/ ao 

a-WoO 















Ao a, 



and 











be the partitions of yl. A, and with respect to the decompositions R x Wo x Wq"*", 
(R X Wo) X W^, and (R x Wq) x W^, respectively. Note that 

and >l(>Co) C Co implies that A^^^owoi'^j ^o) = for all (A, Wo) G R x Wo. Therefore 
Aw±owo = 0- Since A e SO'''(*) we also have * o (A x A) = ^f,, i.e.. 











Ao A, 
A 



Ao A, 
A 



*0 
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This is equivalent to the foUowings. 



(i o (^0 X Ao) = *o- 

(i^) ^oo(^oxA^^x) = 0. 

[Hi) ^'^^x o {Aw± X Aw±) = ^'lyx. 

□ 



3.2. A maximal invariant statistic for testing the hypothesis 

Now consider the action of Go on £ , i.e., the restriction of the (transitive) action of 
G on £ to Go- We define the statistic 



^'•^'•^"^ ^ w„((a:wo),(a.w„)) - 

where (A, Wo, w^) represents a typical element of IR x Wq Wq"*". 

Proposition 3.2. The mapping m in Eq. (9) is a faithful representation of the orbit 
projection tt /Go and therefore a maximal invariant statistic. 

Proof. First for notational convenience set ||wp = ^vi'(w, w), for w e W. Thus 
m((A, Wo, w^)) = — — 2 — V(A, Wo, w^) e £. 

A — II Wo 1 1 

For each {a, A) e Go by partitioning as in Proposition 3.1 we obtain 

||aA^±Wff||2 

m((a,>l)(A,wo,Wo^)) = 



a A) 2 — ||aAvi/oWo| 





I Wo 



±112 



A2- ||wo||2' 

This shows that m is invariant under the action of Go on £. Now suppose that 

m((A, Wo, w^)) = m((/x, uo, u^)), 

i.e., 

Ilw-L||2 llii-L||2 
l|Wo II ||Uo II 



A2-||wo|p /x2-||uo||2' 
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If U|j" = , then Wq — and, since Go acts transitively on £0, we can find {a, A) G Go 
such that 

(a, A) (A, Wo, 0) = (/x,uo,0), 

which shows that (//, Uo, Uq ) and (//, Wo, Wq-) are in the same Go-orbit. 
Now assume 7^ 0, and thereforeUg 7^ 0. We set 



a :- 



1^0 II _ / /"^ - IIuqP 



|w^ll ]j - ||wo||^' 
Choose Aq e 5'0+(*o) and A^^± e S'0(*^x) such that 

(A,Wo) (/^,Uo) 



Ao{ 



- ||Wo||2 - ||Uo||2' 



and 



For A := ( , ) we have (a, A) e Go and 



A^s 



(a, A) (A, Wo, w^) = a(Ao(A,Wo),>lwx(w^)) 



- Wo I . X llw^ll ^ 
-^-^-p(/x,Uo),^-^Uo) 
A* ~ l|Uo|| ||Uq II 

= a(a~^(/i, uq), a~^Uo ) 

= (/i,Uo,Uo). 



□ 



Remark 3.1. From statistical point of view, it is more useful to write the maximal 
invariant statistics m in Eq. (9) as 

™(^.-) = !Sq^ V(A,w)6£. 

where p : W ^ Wq is the orthogonal projection of W onto Wq. 



8 



3.3. Testing the hypothesis (Tl) 

Finally we are in the position to give a test statistic for testing the hypothesis (Tl). 
First note that the hypothesis (Tl) is invariant under Go- 
Theorem 3.1. Consider W^^^^^, the Wishart distribution on L, where cr e £0 o,nd 
i{o) & C is the embedding of a into C Let t : C ^ jCq be the mapping (y,z) i->- 
{y,p{z)), where {y,z) e C Then for the mapping 

{t, m): C £0 X M++ 

(y,z) ^ ((y,p(z)),m(y,z)) 

we have (t,m){W^^^^^^) = (g) m( ) . 

Proof. First write VV^^^^-, as a density with respect to the invariant measure 

du^yjZ) = (y — \\z\\ )' ~d,ydz, 



and rewrite the density in terms of t{y,z) and m{y,z). Then apply [1, Lemma 3] to 
z/ and state the transformation result. □ 

Corollary 3.1. Let mo := dim]K(lVo) (ind mi = m — mo. Then the transformed 
measure tTi(>V^^(^j)) is the beta distribution [5{ri — ^)- 

Proof. We start with writing (i, m) : C ^ Cq x as the composition of the 
mappings 

(y,z) ^ (2/,Zo,Zi) := I l/,p(z), 



^ (y,Zo,ii) := (y,Zo, ||zi|p). 

Using these compositions we transfer the probabihty density of the Wishart distri- 
bution Wf as follows: 

dWi,W,M ^ ^^^(l-|l^iH')'-°*'l(IKII<i)(^i)''<»(!/,z„)<iz, 

r(// — ^) , „ m+1 mi 



r(,7-^)r(^) 

777 

" 2 



d>V,^°(|/,Zo)ci/3(7^-^,^)(77). 
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Thusm(t/,z)~/3(r7-^,^). □ 

Proposition 3.3. The likelihood ratio LR statistic Q for hypothesis (Tl) is given 
by 

Moreover, Q is independent >V^|} and (3{r] - ) . 

Proof. First note that for the submodel (8) and the observation {y,z), the ML 
estimator of a & Cq becomes 

S'i((y,z)) = (y,p(z)), 

where p is the orthogonal projection on the subspace Wq. Thus LR is given by Eq. 
(10). The rest of the proof are direct consequences of Proposition 3.1 and Corollary 
3.1. □ 

Note that the Proposition 3.3 reduces to Theorem 5 in [8] for choice oi rj = N/A. 

4. The Bartlett's test for Lorentz Wishart Models 

As before, suppose t] is known, and consider the statistical model 

and its submodel 

Consider the hypothesis 

Hq : (Tl — (72 — (7 vs. H : ai ^ 02- (T2) 

J^.l. A maximal invariant statistic associated with the hypothesis (T2) 

Note that C can be considered a subcone of £ x £, via diagonal embedding. Also 
the action of G on £ can be, canonically, extended to an action on C x C given by 

g{ai, CT2) = {gai, ga2) ^g e G, V(cti, cts) G £ x jZ. (11) 

Under this consideration the hypothesis (T2) is invariant under G, and therefore a 
maximal invariant statistic is desired. To obtain a maximal invariant we proceed 
with the following lemma. 
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Lemma 4.1. For every cxi and (J2 G £ 



*((Ji,(J2) > V'^(^i^^i)*(^2y^, (12) 
and the equality holds if and only if gci — (T2 for some g & G. 
Proof. Let Ui • (72 := ^(cri,(J2), where cti = (A, w) and (T2 = (/^jU). Then 

(f7i-(72)'-||(7i|n|(72||' = (A/X-W.u)2-(A2-||w|n(/x2-||u||2) 

> A^||u||^ - 2A/i||u||||w|| + n'^Wwf 
= (A||u||-H|w||r>0. 

The equahty holds if and only if A||u|| = //.||w|| which is equivalent to ||c"2|| = t- 11*^1 II- 

A 

□ 

Notation 1. In the remainder of this paper for brevity we use the notation a\ ■ a-i 
for the quadratic product ^'(cri, (T2). 

Proposition 4.1. The mapping n : JC x C ^ ^++ ^^^^ (ci,cr2) ^ (^1)^2); where 



(Ti • (72 + s/iai ■ (72)2 - ||(Ti 



2 



(^1 := — — — and 



(Ti • (T2 - a/ (cTi ■ a2)2 - ||(7i|P||(72 



II 119 ' 

is maximal invariant under the action of G, defined in Eq. (11), and ^1, ^2 o-f^ 
eigenvalues of a-i with respect to a\ . 

Proof. By Lemma 4.1 the mapping tt is well-defined. Next we show that ^1 > ^2 are 
eigenvalues of a\ with respect to Ui. Note that the characteristic polynomial of a-i 
with respect to a\ is 

V{t) ■= det((72 - £ai) = \\a^fe - 2((Ti • a<,Y - \\a<,f . (13) 

One then can easily check that ^1 and ,^2 are indeed the roots of p(^). Next we 
show that TT is, moreover, onto. Let the real numbers ri > r2 > be given. Set 

A := ^ Choose a vector w G such that ||w|| = 2 ^^'^ ^ ~ ^^'^ 
e = (1,0). Then 7r(((7, e)) = (ri,r2) as desired. It is clear that vr is invariant under 
the action of G. Suppose 7r((Ti,(T2) = 7i{a[,a2)- Since G acts transitively on C and 
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TT is invariant under the action of G, without loss of generahty, we may assume that 
(72 = = e. We have 

7r((Ti,e) = 7r{a[,e). 

Eq.(13) imphes that the characteristic polynomials := ||(7i|p£^ — 2((7i.e)£— 1 and 
P2(^) lkilP-^^~2((7i.e)£— 1 are identical.Therefore, if (7i = (A,w) and a[ — (A',u), 
then 

A = ai ■ e — a'l ■ e — X' 

\2 II ||2 II ||2 II/II2 \/2 II ||2 

A — ||w|| = ||(7i|| = ||(T II = A — ||u|| . 
Thus (cTi, e) and (cr^, e) are in the same G-orbit. □ 

4-2. Testing the hypothesis (T2) 

Now by using the maximal invariant obtained in Proposition 4.1 we can test the 
hypothesis (T2) as follows. 

Theorem 4.1. For the observation (ti,T2) & C x C, the ML estimator of a under 
Hq is 

X ^1 + '^2 

<y{ri,T2) := — ^ — , 
and the LR statistic for testing hypothesis (T2) is 

\ 



where > ^2 are eigenvalues o/r2 wi^/i respect to ti. Furthermore, under the hypoth- 
esis Hq the statistics (x(ri,r2) and 7r(ri,r2) = {^1,^2) f^re independently distributed, 
5'(Ti,r2) ~ ^T^CT '^^'^ density o/7r(ri,r2) = (^1)^2) 5'«t'en 6y 

(2^)"-^ ,.2 (n.LiC.o-^ 



Sar?,r^)r^(n-2)^^^ (n,=i(l + ^.O)'^ ' 

Proof. This follows from Proposition 4.1 and Theorem 6.1 in [3] when the irreducible 
cone Q is the Lorentz cone C. □ 

We should mention that the proof of Theorem 6.1 [3], which simultaneously ap- 
plies to all five types of irreducible cones, heavily rests on the analysis of simple 
Euchdean Jordan algebras. 
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Remark 4.1. Recall that in the classical multivariate statistics, the Bartlett's test is 
testing 

Ho-.^^a^In vs. H:E^a''ln, (14) 

for a Gaussian model, where the sample space is R", the distributions are multivariate 

normal distribution AfniO,!^) and the parameter space is PZ}.„(R). Therefore the 
Bartlett's test is testing whether n univariate Gaussian distributions are independent 
and have the same variance a. 

5. Closing Remarks 

In this paper we have shown how maximal invariant statistics can be derived 
and used for testing two specific invariant statistical hypotheses for Lorentz Wishart 
models. Analogs of such hypotheses have been already studied for real, complex and 
quaternion type Wishart models [1]. An interesting topic of future research, which 
its analog has been studied in [1], is testing the hypothesis that the scale parameter 
S of the real Wishart distribution has a Lorentz structure. To clarify what is meant 
by a Lorentz type we note that every Lorenz cone C is isomorphic to a subcone 
of P£)„(R) for a suitable n (see [6], [8] for detail), which means that there exist a 
linear injection p : W x W ^ W^^"- such that p(A,w) is positive definite, for each 
(A,w) G C. Therefore testing whether E has a Lorentz structure requires to show 
that S = p(A, w) for some (A, w) e C. 
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